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Under quite general conditions, we prove that for classical 
many-body lattice Hamiltonians in one dimension (ID) total 
momentum conservation implies anomalous conductivity in 
the sense of the divergence of the Kubo expression for the 
coefficient of thermal conductivity, n. Our results provide 
rigorous confirmation and explanation of many of the existing 
"surprising" numerical studies of anomalous conductivity in 
ID classical lattices, including the celebrated Fermi-Pasta- 
Ulam problem. 
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Since the pioneering work of Fermi, Pasta, and Ulam 
(FPU) revealed the "remarkable little discovery" Q that 
even in strongly nonlinear one-dimensional (ID) classi- 
cal lattices recurrences of the initial state prevented the 
equipartition of energy and consequent thermalization, 
the related issues of thermalization, transport, and heat 
conduction in ID lattices have been sources of continu- 
ing interest (and frustration!) for several generations of 
physicists. The complex of questions following from the 
FPU study involves the interrelations among equiparti- 
tion of energy (is there equipartition? in which modes?), 
local thermal equilibrium (does the system reach a well- 
defined temperature locally? if so, what is it?), and trans- 
port of energy/heat (does the system obey Fourier's heat 
law? If not, what is the nature of the abnormal trans- 
port?) In sorting through these questions, it is important 
to recall that the study of heat conduction (Fourier's heat 
law) is the search for a non-equilibrium steady state in 
which heat flows across the system, but the situation is 
usually analyzed, using the Green-Kubo formalism of lin- 
ear response Q , in terms of the correlation functions in 
the thermal equilibrium (grand canonical) state. A series 
of reviews spread over nearly two decades has provided 
snapshots of the understanding (and confusion) at differ- 
ent stages of this odyssey ||-[| . 

Much of the past effort has been devoted to attempts 
to verify Fourier's law of heat conduction 



(J) = kVT, 



(1) 



where in ID the gradient is replaced by the derivative 
with respect to x. Here, K is the transport coefficient of 
thermal conductivity. Strictly speaking, k is well defined 
only for a system that obeys Fourier's law and where a 
linear temperature gradient is established (for small en- 
ergy gradients such that relative temperature variation 
across the chain is small; in general n is a function of 



temperature, of course). In the literature the dependence 
of k(L) on the size L of the system/chain has also been 
used to characterize the (degree of) anomalous transport. 
However, the definition of n for an anomalous conductor, 
where no internal temperature gradient may be estab- 
lished, is ambiguous. Typically, one defines it in the 
"global" sense, as k(L) = kq = JL/AT, where AT is 
the total temperature difference between the two ther- 
mal baths. However, if the temperature gradient is not 
constant across the system, and/or if there are finite tem- 
perature gaps between the thermal baths and the edges 
of the system due to system-bath contact, one should de- 
fine and study a local k, k = kl = -^f , where VT is the 
local thermal gradient. A very wide range of results have 
been produced by previous studies of different systems: 

• in acoustic harmonic chains, rigorous results es- 
tablish that no thermal gradient can be formed in 
the system, with the result that formally kq ~ L , 
which can be understood heuristically by the sta- 
bility of the linear Fourier modes and the absence 
of mode-mode coupling; 

• in the "Toda lattice," an integrable lattice model 
[p|JlO|, in which the result kq ~ L 1 (llj can be 

understood in terms of stable, uncoupled nonlinear 
modes, the solitons, which are a consequence of the 
system's complete integrability ; 

• in non-integrable models with smooth potentials, 
including (i) the FPU system, leading eventually 
to claim that chaos was necessary and sufficient for 
normal conductivity (kq — kl ~ L°) a claim 
that has been countered by convincing numerical 
evidence for anomalous conductivity in FPU chains 
(kl ~ L 4 ) [jl2|,^3); (ii) the diatomic (and hence 
non-integrable) Toda lattice, where initial numer- 
ical results claiming kl ~ L° [[id} have recently 
been refuted by a more systematic study showing 
n L ~ L 0A Q; and (iii) the "Frenkel-Kontorova 
model," where recent studies have shown that (at 
least for low temperatures) kl ~ L° ; 

• in non-integrable models with hard-core potentials, 
including (i) the "ding-a-ling" model p7[ ]; (ii) the 
"ding-dong" model [Q, and (iii) even simpler sin- 
gle particle chaotic billiard model 19 1, where nu- 
merics show convincingly that kl — kg ~ L ■ 

This bewildering array of results has recently been par- 
tially clarified in a scries of independent but overlapping 
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studies. The numerical studies of Hu et al. |T(| and of 
Hatano fl5|| show that overall momentum conservation 
appears to a key factor in anomalous transport in ID 
lattices. Lepri et al. and Hatano |TJ|] have ar- 

gued that the anomalous transport in momentum con- 
serving systems can be understood in terms of low fre- 
quency, long-wavelength "hydrodynamic modes" that ex- 
ist in typical momentum conserving systems and that 
hydrodynamic arguments may explain the exponents ob- 
served in FPU |t],|l| and diatomic Toda lattice |ij§. 

In the present Letter, we extend and formalize these 
recent results and resolve finally at least one important 
aspect of conductivity in ID lattices: namely, we present 
a rigorous proof that in ID conservation of total momen- 
tum implies anomalous conductivity provided only that 
the average pressure is non-vanishing in thermodynamic 
limit. 

We consider the general class of classical ID many- 
body Hamiltonians 



N-l , 1 v 

H =^2 [ n Pn + K+l/2(?n+l - qn)) , 

n=0 \ Zl7ln > 



(2) 



where V n+ 1/2(9) is an arbitrary (generally non-linear) in- 
terparticle interaction. Note that the potential, T^ n +i/2, 
depends only on the differences between two adjacent 
sites; in particular, there is no "on-site" potential, 
Uos(qn), depending on the individual coordinates. The 
(finite) system is considered to be defined on a system 
of length L — Na with periodic boundary conditions 
(<7jv iPn) = (<Zo,Po), where the actual particle positions 
are x n = na + q n . In our analysis the masses m„, as well 
as interparticle potentials V n+ i/ 2 (q), can have arbitrary 
dependence on the sites n, though the examples stud- 
ied in literature to date have mostly had uniform poten- 
tials V n+1 / 2 (q) = V(q) and uniform, m n — m, or dimer- 
ized m-m = rai,m2 n +i — ^2, masses. We require only 
that the Hamiltonian (^) be invariant under translations 
q-n —* q-n + b for arbitrary b. This requires J7os(<7n) = 
Ji7j[ . Note that we may write the Hamiltonian in Eq. (|^) 



as H = En=o h 
density 



-1/2 



l n+l/2 



2 

Pn+l 



where h r , 



+ -r^ + V n+1/2 (q n+1 -q n ), (3) 



-1/2 



is the Hamiltonian 



4m n+ i 4m n 

Our aim is to estimate n, the coefficient of thermal 
conductivity, which is given by the Kubo formula [E2J 



lim 



lira — 

L — >oc L 



dt(J{t)J} . 



(4) 



Here we have written the canonical average of an 
observable A at inverse temperature (3 as (A)p = 
/ Il n dp n dq n Aexp(-f3H)/J Tl n dp n dq„exp(-f3H). The 
order of limits in Eq. ([l|) is crucial to the precise def- 
inition of K [§. In Eq.(|4|), J = J2 n =o 3 n 1S tne total heat 
current, and j n is the heat current density Eq|, given by 



jn — {h n +l/2, hn-l/2} — (5) 

= - q n ) + V n _ 1/2 (q n - g„_i)J , 

where {, } is the usual canonical Poisson bracket. 

Using Eq.s (|^,||), we find that current density given by 
(o) satisfies the continuity equation 



^n+1/2 — {H, Ki+1/2} — 3n+i - j 



(6) 



Our ensuing analysis is similar to that used by Mazur 
p3|, with a crucial difference: we will average correlation 
functions over a finite rather than infinite time domain, 
T. We start with an elementary inequality. For an arbi- 
trary observable X(t) = X ({q n (t) , p n (t)}) , we have 



dtgr(t)(X(t)X)p > 



(7) 



where gr(t) is a suitable L 2 (M.) window function of effec- 
tive width T, which has the following properties: 

(i) J™ oo dtg T (t)=T. 

(ii) JZ dtg 2 T (t)=T. 

(hi) g(u) := dtg T {t)e wt > for all w. 

The natural choice satisfying these conditions is a Gaus- 
sian, gr(t) = v2 cxp(— 2it{t/T) 2 ). Using elementary 
Fourier analysis, the above inequality (f7j) is easily proved 
by rewriting it as 



dujg T (uj)(Sx{uj))p > 



(8) 



where Sx(tu) = liniT- 



i_ 

->oo y 



fidte™' 



X(t) 



is the power 

spectrum of the signal X(t). Obviously, Sx(u>) > 0, and 
given (iii), the inequality (0J|) is clearly fulfilled. We 
now write the observable X as X = A + aB, a € K. 
Optimizing with respect to the parameter a, we arrive 
at the Schwartz-like inequality 



dtg T {t){A{t)A) p / dtg T (t)(B(t)B) p > 



> / dtg T {t){B{t)A) p 



(9) 



The above inequality is of quite general use. We im- 
plement it by taking A = J and B = P, where P = 
12n=o Pn is the total momentum. For Hamiltonians of 
the form (||), P is an integral of motion P = {H, P} = 
due to translational symmetry. Since P(t) — P, the in- 
equality (0) reads 



dtg T {t){J(t)J)p > T 



(JP)l 



(10) 



The RHS of Eq.^) can be easily evaluated: (P 2 )p = 
N/(3, and (JP) P = /T 1 ^"=0 (V (q n +i - Qn))p 
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since we have in general that (A({q n })B({p n }))p = 
(A({q n }))p(B({q n }))f3. (V^ +1/2 (q n +i - q n ))p is an aver- 
age force between particles n and n + 1, i.e. the thermo- 
dynamic pressure and does not depend on n. (In thermal 
equilibrium the average net force on particle n vanishes 
and hence (V^_ 1/2 (q n -q n -i))/3 = (U,' +1/2 (<?n+i -q n ))p-) 
The pressure can be rewritten through the usual thermo- 
dynamic definition 

± _dF 1 OF 1^',,, 

^=~dL = N~da~ = L^ \ V n+i/2{ x n+i -x n + a))p, 

where exp(— (3F) = J U n dp n dq n exp(—(3H). Inserting 
the above and multiplying with {3/L, we find that in- 
equality (||) reads 

| J dtg T (t)(J(t)J) f3 > T(/) 2 . (11) 

Since the Kubo formula can be equivalently written in 
terms of the win- 

dow function as n = 2~ 1 / 2 limT^ 00 J dtgT(t)C(t), where 
C(t) = lim i _ >00 (J(t)J) | g/L, since lim^oo g T (t) = y/2, 
and implementing in the above result (|ll]) the two limits 
as indicated in (|j), we have proved our main result: 

Theorem: In momentum conserving systems of type 
, if the pressure is non-vanishing in the thermody- 
namic limit, limi^oo <j) > 0, then the thermal conductiv- 
ity diverges and n — > oo. 

Therefore, we find anomalous energy transport as a sim- 
ple consequence of the total momentum conservation. 
The only case in which the pressure is expected to vanish 
at any temperature is when the forces between particles 
at zero temperature equilibrium are zero (L^ +1 / 2 (0) = 0) 
and the interparticle potentials are all even functions 
(V n+ i/ 2 (q) = V n+ i/2 (—9)) so that the forces are also ex- 
pected (and found numerically for (3 FPU problem) to 
average to zero for arbitrary canonical thermal fluctua- 
tions. This is indeed the case for the (3 FPU problem, 
where V n+1/2 (q) = \q 2 + \f3q A and there the 

integrated correlation function diverges for more subtle 
(dynamical) reasons (the slow asymptotic power-law de- 
cay of current-current correlation function ~ t -0 6 ). 

Even if the zero temperature equilibrium forces vanish 
V^'+i/2(0) = 0, we still have non- vanishing finite temper- 
ature pressure (due to 'thermal expansion' of a system 
confined to a fixed volume L = aN) whenever interpar- 
ticle potentials are not even. This is the case for the a 
FPU model, V n+ i/ 2 (q) = \q 2 + \aq 3 , for the modified 
diatomic Toda lattice [pi, V n+ i/ 2 (q) = cxp(— q) +q, and 
for the diatomic hard-point ID gas |24|,[l5| V n+ i/ 2 (q) = 
{0 if q > —a; 00 if q < —a}. For the usual diatomic 
Toda lattice [p|, V n+1 / 2 = cxp(— q) the pressure is non- 
vanishing even at zero temperature, since V^ +1 / 2 (0) 7^ 0. 



To augment and illustrate our analytic discussion, we 
have simulated numerically the current-current autocor- 
relation function (J(t)J)p/L in a generic anharmonic 
chain, namely in the "a/?" FPU model with V n+ i/ 2 (q) = 
\q 2 + ^aq 3 + j(3q 4 where we take a = 2, (3 — 4 and 
energy per particle E/N = 1. In Fig. 1 we compare 
the results for N — 16, 32, 64 with the equillibrium value 
of the squared pressure 2 = 0.964 .... We have also 
checked numerically that for the symmetric interparticle 
potential (same as above except with a = 0) the pres- 
sure indeed vanishes and the current-current correlation 
functions decay asymptotically as ~ t~ 6 (in agreement 
with results of Refs. 

Given that momentum conservation implies anomalous 
conductivity, it is natural to ask whether the converse is 
true: namely, does anomalous conductivity imply that 
the model conserves momentum? Two counterexamples 
show that this result is not true. First, if one considers a 
linear chain of optical phonons — soU„ +1 / 2 ~ (?n+i — 9n) 2 
and Uos ~ In — one can show E3j by a straightforward 
extension of the arguments of Ref. that this mo- 
mentum non-conserving model nonetheless has anoma- 
lous transport. Similarly, there is a momentum non- 
conserving integrable model due to Izergin and Korepin 
p6| that also shows anomalous conductivity J25| . Finally, 
let us to stress that in ID lattices the nature of dynamics, 
whether it be completely integrable, completely chaotic, 
or mixed, does not affect our result: if total momentum 
is conserved and the canonical average of the pressure 
does not vanish, the transport is anomalous. We shall 
address the central issue of the necessary and sufficient 
conditions for normal transport in a forthcoming paper 
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FIG. 1. Current-current autocorrelation function in the 
"a/?" FPU model with a = 2, /3 = 4, and E/N = 1. We show 
numerical data canonically averaged over 500 pseudo-random 
initial conditions for three different sizes N = 16, 32, 64 and 
compare it to the squared pressure (f) 2 (dotted line). 
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